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1 Introduction 



In this article we describe a gauge invariant way to compute the one loop 
corrections to processes involving the production and decay of unstable ele- 
mentary particles, such as the W, Z and Higgs bosons or the top quark. All 
these particles have widths which are (possibly) a sizable fraction of their mass 
(r/m ~ 1/40 for the W and Z bosons). A fixed-order amplitude in perturbation 
theory respects the symmetries of the theory, like gauge invariance. However, 
in a reaction which involves the production and decay of these unstable par- 
ticle this amplitude contains a non-integrable infinity due to the propagator 
l/(p 2 — m2 )- We are thus forced to take into account higher order diagrams: 
the (one-particle reducible) self energy graphs of the unstable particle. Resum- 
ming these gives rise to a pole which now is off the real axis, and the amplitude 
is finite. However, one can no longer assume that this amplitude is still gauge 
invariant, as one has mixed different orders in perturbation theory. Experience 
has shown that gauge dependent terms can sometimes be much larger than the 
physical answer, so it is necessary to investigate this problem. 

There are many ways in which the gauge invariance of the amplitude can be 
restored. The easiest way to circumvent the problem is to use the narrow width 
approximation, that is, to factorize the process into production and decay of 
the unstable particle. This is often not good enough when the width of the 
unstable particle is sizable. Another approach is to introduce the width in the 
Lagrangian, by adding and subtracting suitable terms, thus using a complex 
width everywhere. This is a problematic procedure, as we will discuss in section 
[2]. Finally, it is often possible to obtain a gauge invariant amplitude by analyzing 
a particular reaction in detail; this is precisely what has frequently been done in 
Z physics at LEP I. In this paper we present a more general method to obtain a 
gauge invariant amplitude based on the analytical properties of the amplitude. 

This procedure to treat unstable particles has been given a long time ago 
by M. Veltman [?]. Recently it has been applied to the production and decay 
of a single uncharged particle, the Z boson, to first order, in Refs [?, ?, ?]. 
The method centers on the isolation of the position and residue of the unstable 
particle pole from a fixed-order calculation. We generalize it to the treatment 
of multiple, charged unstable particles. Both extensions bring new effects at the 
one loop level. A side effect worth mentioning is that the non-resonant terms, 
which are not enhanced by the resonant propagator, are separated in a gauge 
invariant way. This means that they can be computed to one order less than 
the resonant terms. This greatly economizes on the amount of work necessary. 
The reaction we have in mind is of course W pair production at LEP II; in this 
paper we present the general formalism. 

The layout of this article is as follows. First we review some alternative 
methods which have been proposed to treat unstable particles in perturbation 
theory. Next we separate the diagrams for a given reaction into three groups: 
factorizable resonant, non-factorizable resonant and non-resonant. The first 
group is analyzed to all orders for a single uncharged particle in section |j. A 
list of subtleties is discussed in section g this includes the effects of a charged 
particle. Next we analyze the non-factorizable diagrams. Finally we recapit- 
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ulate the gauge invariance arguments, and give a recipe for tree level and one 
loop calculations in this scheme. The appendices contain various proofs needed 
in the text. 

2 The Finite Width and Kinematics 

The consistent treatment of unstable particles cannot be implemented on the 
level of the Lagrangian. A proposal using this idea [?] is to add and subtract 
a piece proportional to the width in the mass term of the original Lagrangian: 
m 2 cj) 2 — ► (m 2 — imT)(j) 2 + imTcft 2 . The first term now gives a complex mass 
M 2 = m 2 — imT in the propagator, the second one is taken as a perturbatively 
treated two-particle interaction. In the standard model this leads to a complex 
weak mixing angle sin 2 9\y = 1 — Myy/Mg, as MjL and M§ are both complex. 

The major drawback is that this Lagrangian gives rise to unwanted effects 
in case the unstable particles are not near their mass shell. One example would 
be deep inelastic scattering, where the t-channel W boson obtains an on-shell 
width Tw, in spite of the fact that the self energy is purely real for spacelike 
p 2 . Another example is the Bjorken mechanism for Higgs production: e + e~ — ► 
Hfi + [i~ (see Fig. |]). The two diagrams may look similar, but have to be treated 
differently. At ^/s « mz (LEP I) it is necessary to resum the first Z boson, but 
not the second one; at y/s > m# + mz the second one can be on-shell. Using 
the on-shell width in the first Z boson propagator is in fact wrong: as the self 
energy scales roughly with s one should use a width which is four times larger 
than the on-shell width at LEP II energies. This error can be corrected in this 
scheme, but this requires a one loop calculation (introducing errors of C(a 2 )). 




Figure 1: The Bjorken process at y^s rj mz (left) and \fs > mz + (left). 
Z* denotes an off-shell Z boson; the thick line indicates a propagator which 
needs to be resummed. 

The treatment of the unstable particle propagator depends therefore criti- 
cally on the kinematical configuration. We will refer to a particle as unstable 
only when it can kinematically be on its mass shell. In this case we are forced 
to resum its propagator. Otherwise the particle can and should be treated 
perturbatively. 

The basic idea of the method used here is to consider the analytical proper- 
ties of the transition amplitude. Neglecting for the moment all further details 



2 



and qualifiers, the amplitude containing one unstable particle has the form 

in which p is the momentum flowing through the unstable particle propaga- 
tor, M 2 the (complex) pole position, w the residue at this pole and n{p 2 ) the 
non-resonant remainder. The pole position, the residue at the pole and the 
remainder are gauge invariant quantities. We will show how these can be com- 
puted systematically and show this gauge invariance explicitly. Initially we will 
discuss this for a single particle. The generalization to more unstable particles 
is given afterwards. 

3 Classes of Diagrams 

In order to simplify the discussion we will divide the diagrams — of a per- 
turbative expansion of the amplitude to all orders — into three classes. This 
separation is not gauge invariant; we use it to extract the gauge invariant con- 
tributions to the residue at the pole in Eq. (|l]) more easily. 




production propagator decay 
Figure 2: The structure of the factorizable resonant diagrams 

The first class of diagrams we call the factorizable resonant diagrams (see 
Fig. ||). It comprises all diagrams that retain the original unstable particle prop- 
agator l/(p 2 — m 2 ) outside the loops. These are the diagrams appearing in the 
narrow width approximation. The radiative corrections occur separately in the 
production, the decay and the propagator, but do not connect production and 
decay (or, in the case of more resonances, the decays of different resonances). 
The contribution to the amplitude can be written as 

"4fact = o T~ I — 2 2 ) ' (2) 

where m denotes the renormalized mass, which is real and finite, but need not 
be the on-shell mass (which is defined as the real part of the pole position). The 
corrections to production and decay are contained in W, and II (p 2 ) is the one 
particle irreducible self energy. This part of the amplitude is infrared divergent. 
The structure of the soft bremsstrahlung diagrams which have to be added to 
cancel this divergence is identical (after integrating out the extra soft photons). 

The ellipsis in the argument of W stands for other variables which deter- 
mine the kinematical configuration and hence the the amplitude. Examples are 
the angles of the resonance decay products and the momentum squared of par- 
ticles produced far off their mass shell (for instance the second Z boson in the 
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Bjorken mechanism for Higgs production at ^/s ~ mz)- It is important that 
the integration boundaries of these variables are independent of p 2 , unlike, e.g., 
Mandelstam invariants. This would introduce extra terms in the expansions 
used below. On the other hand, angles in the CMS are obviously a good choice. 
If this does not suffice to eliminate all variables with p 2 dependent integration 
boundaries (like the momenta squared mentioned before) one has to introduce 
explicit mappings from variables < X{ < 1 to these p 2 dependent variables. 
The angles and mappings are not unique; one can for instance change the frame 
in which the angles are defined. Different choices will change the answer in a 
fixed order calculation, but the difference will be of higher order. An explicit 
proof for the Oip) case can be found in appendix ||. 

This dependence on other variables introduces problems when the first ar- 
gument of W is such that the kinematical configuration becomes unphysical. 
This happens when p 2 is below the threshold for production or decay of the un- 
stable particle. We choose to define W(p 2 ) to be zero in this region. However, 
this will cause problems near, and especially below threshold. 



Figure 3: The structure of the non-factorizable resonant diagrams 

There is a second class of diagrams which also diverges linearly for p 2 — » m 2 . 
This group consists of all infrared divergent (photonic and gluonic) corrections 
to the lowest order process which span production and decay (or the decay of 
different unstable particles). A prototype for this kind of diagram is given in 
Fig. ||. The linear divergence for p 2 — ► m 2 is related to the infrared divergence 
in that both occur in the limit that the photon momentum in the loop Q — > 0. 
However, whereas the infrared divergence is logarithmic, the on-shell divergence 
is linear (or higher for more unstable particles). In both cases the divergence 
cancels against the corresponding soft Bremsstrahlung graphs (which have to 
be treated analogously to the virtual graphs), but only after integrating in the 
soft-photon approximation the photon up to an energy much larger than the 
width. 

Thirdly, all other diagrams are non-resonant, i.e., they do not diverge in 
the limit p 2 — ► m 2 . Examples of these are diagrams leading to the same final 
state without the resonant propagator, and all non-infrared divergent graphs 
connecting production and decay (or different decays of unstable particles). 
These graphs can not and should not be resummed. 

In this classification scheme we assume that there are no two channels involv- 
ing different resonances leading to the same final state. An example violating 
this is e + e~ — ► [x^ ix~vv, which can proceed via two W bosons or two Z bosons. 
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These two contributions peak in different regions of phase space, but there is 
an overlap. We do not address this issue in this paper. 



4 The Factorizable Resonant Diagrams 

In this section we discuss how to find the position of the pole and the contri- 
bution to its residue from the factorizable diagrams. Both these quantities are 
gauge invariant. Outside the resonance region p 2 ~ m 2 , Eq. @ can just be 
expanded perturbatively. However, near the pole p 2 = m 2 a fixed-order pertur- 
bative calculation is useless; instead we will have to resum part of the higher 
order self energy contributions. This introduces problems in the case of charged 
unstable particles. In this section we assume that the resonance is neutralfj; the 
extension to charged particles will be given in the next section. 

After summing the geometrical series we isolate the pole structure as [?, ?] 



W(p 2 



p2 _ m 2 _ n(p 2 



W(p 2 ) W(M 2 ) 1 



p 2 - m 2 - n(p 2 ) p 2 — m 2 i — n'(M 2 ) 

w 

p 2 - M 2 



W(M 2 



p 2 — M 2 1 — U'{M 2 ) 



= n(p>) + -^, (3) 



with the pole position defined by M 2 — m 2 — II(M 2 ) = 0. (If m 2 is the on-shell 
mass M 2 = m 2 — imT. This differs from the standard on-shell renormalization 
scheme by the fact that we use II(M 2 ) in the counterterm, rather than II(m 2 ); 
the difference is of C(a 2 ).) Here, w and n(M 2 ) are gauge invariant. However, 
the quantity W(M 2 , . . .) is not yet well-defined because of the dependence on 
the other arguments: if p 2 is off the real axis the angles also have to be taken 
complex. However, we do not need to evaluate W(M 2 , . . .), as the residue can 
be read off from an alternative expansion of the first class of diagrams around 
the real mass m 2 . This gives 

pZ _ m z d—d^ ^pZ _ m z j pZ _ m z z-^ i^pZ _ m zyi 

The coefficients can be found as a series in a ~ II by a straightforward expan- 
sion; the first two terms are given by 

W-i = W(m 2 ) + ^ 2 [W{p 2 )Tl(p 2 )\ p2=m2 + \^^ 
W{p 2 ) - W{m 2 ) 

p2 _ m 2 

W(p 2 )U(p 2 ) - W(m 2 )U(m 2 ) - (p 2 - m 2 )^[W(p 2 )U( P 2 )] p2=m2 
+ ,„ 2 ^2 + < 6 ') 



N(p 2 



(p 2 — m 2 )' 



(In the last term we have chosen not to continue the expansion in positive 
powers of (p 2 — m 2 ), obtaining differences rather than derivatives.) As was 

1 that is, neutral with respect to any unbroken gauge group. 
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shown in Refs [?, ?], the two expressions (H) and (Q) correspond rather nicely 
in low orders in a and (p 2 — to 2 ). In fact, it turns out that 

w = W^i (7) 
n(p 2 ) = N(p 2 ) (8) 

to all orders, and thus the residue at the pole l/(p 2 — M 2 ) and the non-resonant 
remainder can be obtained directly from a perturbative off-shell calculation. A 
proof is given in appendix [A|. To the best of our knowledge, these relations have 
not been stated in full generality before. They permit us to find the residue at 
the pole and the non-resonant parts of the resummed expression ([|) from the 
unresummed one (||). The latter is a straightforward perturbation expansion. 

Let us discuss the terms needed for a one loop calculation and their inter- 
pretation. The first order result is given by 

■4'°' = 5^ W 
p 2 - M 2 w 

M 2 = TO 2 + n (1) (M 2 ), (10) 

integrated over the full phase space. We have dropped the subscript factor- 
izable, as the non-factorizable and non-resonant terms do not give a contri- 
bution in this order. The superscripts indicate to which order the different 
parts are computed: we need the one loop self energy at the pole (in fact 
only the imaginary part) and the on-shell amplitude at tree level. The limit 
r = iH(M 2 )/m — > in the propagator and phase space as well gives the narrow 
width approximation. 

In the next approximation one has to use the two loop corrected on-shell self 
energy Il( 1+2 )(M 2 ) (which gives the physical width) in the definition of the pole 
position in the previous term. Besides this one finds the following contributions. 

- The 0(r) corrections 

Afr> _ jv«V) = wm ^\ - T ( "' 2) ■ (ID 

These corrections (and the non-resonant contributions) have been studied 
in Ref. [?]. 

- The one loop on-shell correction with resonant propagator 



(1) _ w-'{m 
^fact ~ p 2 _ m 2 ' 1 ' 
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- The derivative term 



A ~ p 2_ M 2 • W 

In the limit r — > this term compensates for the fact that, in the nar- 
row width approximation, one has twice a field renormalization of the 
resonant field (Z^) 1 ! 2 = — IT(to 2 ). It also implements in a pertur- 

bative way the p 2 -dependence of the self energy, which is included in most 
descriptions of the Z pole. 
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- The term 

A(W >) _ W^'{m^){w?) 

A ~ p>-M 2 ■ [ > 

(which vanishes in the limit r —* 0) compensates at this order in a the 
dependence on the choice of the other kinematical variables mentioned in 
section |2| A proof is given in appendix [b]. 

In the case that there are more unstable particles which can simultane- 
ously be on their mass shell, the above resummation can be performed for each 
pole separately. The resulting series provides a natural way to isolate the rele- 
vant contributions. In the case of two unstable particles these are the doubly 
resonant terms (both unstable particles on their mass shell) for a tree level cal- 
culation, and the one loop doubly resonant and tree level singly resonant terms 
for the 0{oi) corrections. The non-resonant terms only occur at 0(ct 2 ) (with 
the possibility of large logs). 



5 The Fine Print 

The formulae in section |3] look concise, but several subtleties have not been 
addressed; we will treat them in this section. We first investigate the thresh- 
old behavior and show that the method developed here breaks down in this 
region. Next we consider the definitions of W and II in greater detail, and 
treat the phenomena of mixing and unphysical particles. Finally we expose the 
additional problems that occur when there are charged unstable particles. The 
difficulties with the proper definition of the kinematical variables have already 
been addressed in section |2| and appendix ||. 



5.1 Threshold problems 

The expansions @ and (Q) introduce difficulties near the threshold for the 
production or decay of the unstable particle. In the case of the Z peak this is 
not a problem, as the CMS energy squared s coincides with p 2 . However, this is 
not the case when more particles are produced, for instance e + e~ — ► Hfi + fj,~ at 
s ~ (mf/+mz) 2 via the Bjorken mechanism. The leading term in the expansion 
now is W(m'z, . . .); but if s < (mu + mz) 2 , the point p 2 = m 2 z is outside the 
physical region. The angles between the Higgs and Z bosons would have to 
be complex to satisfy energy-momentum conservation. As stated before, we 
choose to define W(p 2 , . . .) = for unphysical p 2 to avoid these problems. This, 
however, effectively annuls the whole resummation scheme and one finds back 
the original l/(p 2 — m 2 ) divergence in N(p 2 ) as the threshold is approached from 
below. Far below threshold no resummation is needed and no problems occur. 
The bad behavior will also manifest itself when approaching the threshold from 
above due to the occurrence of terms T/ (distance from threshold). Therefore, 
with the choice W(p 2 ) = for unphysical p 2 one cannot use the expansion very 
close to threshold. One should thus either not use the separation in resonant 
and non-resonant contributions (which is possible for the Bjorken process and 
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has been done until now [?, ?, ?]), or attempt an expansion in other variables 
or an non-relativistic approach [?]. 

5.2 The definition of the self energy II 

The split between the self energies II (j? 2 ), which are resummed in Eq. (||), and 
the corrections in W(p 2 ) to production or decay, is not unique. Pieces of II(p 2 ) 
that vanish at the pole position, and thus do not contribute to the physical 
width, can also be assigned to the perturbatively treated part W(p 2 ). In fact, 
one can choose to treat all p 2 dependence in the self energy perturbatively 
which implies adding terms proportional to II(p 2 ) — II(M 2 ) to W(p 2 ). In this 
way one can achieve H'(m 2 ) = and the equivalence of the two series (^) and 
(||) is trivial. 

In the case of the Z and W bosons it is customary to resum the p 2 -dependent 
fermionic self energy. This is a gauge invariant quantity, even off-shell. The 
bosonic self energy does not contribute on-shell and is calculated perturbatively. 
This perturbative treatment is required by the gauge cancellations with vertex 
corrections and boxes. The scheme given here is different, in that we only resum 
the physical width, and also treat the p 2 dependence of the fermionic self energy 
perturbatively. The difference is always of higher order in a. 

5.3 Mixing and unphysical particles 

At one loop and higher there are graphs converting physical particles into each 
other, such as the photon and Z boson (see for instance [?, ?, ?, ?]). The 2x2 
matrix for the self energies can easily be diagonalized; the resulting eigenvectors 
define the n loop photon and Z propagators. The Z propagator can then be 
resummed as before. 

In a general gauge there also is a plethora of unphysical particles included 
in the standard model Lagrangian: Fadeev-Popov ghosts and the unphysical 
Higgs bosons. These have masses which depend on the gauge parameter £, 
and it is not clear how these should be resummed, although various suggestions 
have been made [?]. As we isolate the pole structure from a fixed order off- 
shell calculation with real masses via Eqs ([?]) and (||) all gauge dependence has 
already disappeared before the resummation is performed. Only the physical 
width, which cannot depend on the gauge parameter, is resummed and one 
never encounters these unphysical particles in the parts taken from the higher 
order graphs. 

5.4 Charged particles. 

When the unstable particle is (electrically or color) charged, a complication 
arises from logarithmically divergent on-shell graphs. These are graphs that 
contain terms proportional to log(p 2 — m 2 ) because the massless particle cou- 
ples to the resonance. Examples of one loop graphs of this type are given in 
Fig. f|; these graphs would be infrared divergent if the unstable particle would be 
on its mass shell. This means that WW (m 2 ) is not well-defined (not regulating 
this divergence with a photon mass), invalidating the expansion (Q). One could 
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argue that these on-shell divergent logarithms should cancel against the corre- 
sponding soft bremsstrahlung graphs, as p 2 — m 2 ^ is just another infrared 
regulator, on which the physical result should not depend. Unfortunately this 
only holds in the limit E~ < cu <C all other scales (to validate the soft-photon 
approximation), but lo 3> T (so that V is only an infrared regulator). In case 
the width is non-negligible these two requirements are clearly incompatible. 




Figure 4: The two on-shell divergent scalar three point functions with physical 
examples 

We can still use the expansion techniques by defining the function W^(m 2 ) 
in the representation of the amplitude as a sum of fundamental loop integrals 
times coefficient functions. Fundamental integrals are a set of linearly inde- 
pendent integrals in which all tensor integrals can be expressed. In the case 
of one loop corrections, these are the scalar one loop integrals with non-zero 
kinematical determinants, but for higher loops there are also fundamental ten- 
sor integrals which cannot be algebraically reduced to scalar integrals. A proof 
that in a one loop amplitude the scalar integrals are the fundamental integrals is 
given in appendix [C]. This independence implies that the coefficient functions of 
physical fundamental integrals (integrals that depend only on physical masses) 
are gauge invariant, and the coefficients of unphysical fundamental integrals 
(those containing a gauge parameter £) must vanish. Now, the log(p 2 — m?) 
terms only occur in the scalar functions, as the coefficient function is a ratio- 
nal function. We are therefore free to include finite width effects in the scalar 
functions without concern about gauge invariance. 

An obvious approach is to apply the expansion (||) to the unstable particle 
propagators before the integration over the loop momentum Q when evaluating 
these scalar functions. Up to O(a), the complicated non-resonant terms in Eq. ||| 
do not contribute. In the simplest case, the vertex corrections, one just obtains 
the scalar three point function with complex masses. The evaluation of this 
function was already described in Ref. [?]. A numerically stable implementation 
is given in Refs [?, ?]. Note that one should only use a complex mass in 
propagators which are resonant in the infrared limit Q — > 0, as in the other 
propagators the expansion of the self energy around p 2 = M 2 does not describe 
its behavior well in the dominant region. The extension to higher order graphs 
would be much more involved. 

The higher point functions can now also easily be evaluated. One just 



9 



subtracts from the infrared divergent n point function with real masses the 
corresponding infrared divergent three point function (suitably scaled^), and 
adds the same three point function with complex masses. For the four point 
function this is just the first term of the Taylor expansion proposed in Ref. [?].0 
This way also the large cancellations between three and four point functions 
near the edge of phase space ( | cos 9 e w | w 1 in the case of W pair production) 
are unaffected by the width. Note that on the edge the four point function 
reduces to a sum of three point functions. 

The prescription followed for the scalar integrals must be matched, of course, 
by the prescription used in the soft bremsstrahlung integrals; the lower end of 
which should be the pole at Q 2 = of the scalar functions. 

w 

Figure 5: Example of a diagram which is not analytic at p 2 = m 2 : the one loop 
photonic contribution to the self energy of a W boson. 

There is one more graph which spoils the expansions; this is the self energy 
graph of Fig. [|, which also is non-analytic at p 2 = m 2 . The structure of this 
diagram is given by 



2 2 B (p 2 ,m 2 ,0) - B (m 2 ,m 2 ,0) 
9l{P ) + 92{P ) 2 — — 



with the gi(p 2 ) non-singular functions, Bq the one loop scalar two point function 
and 

B (p 2 ,m 2 , 0) - B (m 2 ,m 2 , 0) 



p 2 — m 2 




(16) 



A resummation of the propagators before integration therefore gives a complex 
mass in the one remaining overall propagator in Eq. (|l^) and in the internal 
propagator, which just means replacing m 2 — > M 2 everywhere in Eq. (^). This 
also shifts the subtraction point from m 2 to M 2 , as required in section |j. 

5.5 Further refinements 

As mentioned above, the details of the evaluation of the scalar functions do 
not influence the gauge invar iance of the final result. For any non-singular 
approximation of the scalar function we have a gauge invariant residue. We 
will mention some variants of the method just described. 

2 One can either multiply by the non-infrared propagators, or the ratio of overall determi- 
nants occurring in the scalar functions. 

3 This article contains an error, in that to obtain a convergent Taylor series it is not sufficient 
to subtract the non-analytic parts, but one has to replace these by the off shell expressions. 
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- One can use complex masses in the divergent logarithms only, rather than 
when resonant in the scalar function. The difference is formally of order 
aTf/j, only, with [i some scale. However, this scale may be related to the 
distance from threshold rather than the mass of the resonance, so in that 
case one may miss some important threshold effects. An example is the 
Coulomb singularity, which would not be smeared by the width. 

- One can use a complex mass everywhere in the scalar functions. As 
mentioned above, this introduces errors, but these are of one order higher 
in a and thus irrelevant. 

- There are terms (p 2 — m 2 ) \og(p 2 — m 2 ) associated with the non- infrared 
radiation of a photon off a charged unstable particle. As these also origi- 
nate at Q ~ their resummation to (p 2 — M 2 ) \og(p 2 — M 2 ) rather than 
putting them to zero should improve the resonant behavior. An example 
of a diagram in which this occurs is given in Fig. ^. 




Figure 6: Example of a subdiagram with a photonic threshold. The (predomi- 
nantly) t-channel W propagators are of course not resummed. 

- Similarly, at the threshold for on-shell production of two particles there 
is a term 

v 7 -^ , ( m 2 + ml -p 2 - 2\f^5 
log 2 — v - 

with 6 = m 2 m 2 — (p 2 — m 2 — m 2 ,) 2 /^ related to the Kallen function. At this 
threshold, the scalar function is also dominated by Q ~ 0, so including 
a complex mass for these particles should give a better approximation. 
This happens for instance in the diagrams of Fig. Q with the 7 replaced 
by a Z and q 2 ~ (mi + m.2) 2 . 

- It is not necessary to put p 2 = m 2 in the non-divergent parts of the 
evaluation of the scalar functions; in fact, leaving the momenta off shell 
(but the masses real) will also improve the threshold behavior. 

6 The non-factorizable resonant diagrams 

We now turn to the non-factorizable resonant diagrams. They are absent at 
tree level. We will discuss them in three steps. The first one is to separate 
the resonant contribution of these diagrams from the non-resonant part. The 
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non-resonant part can easily be evaluated, but the on-shell divergence of the 
resonant terms will have to be resummed. After showing how this can be done 
for one loop integrals we note the explicit structure of the divergent terms 
in the integrals needed for e + e~ — * 4 fermions. Finally, we comment on the 
cancellations between these diagrams and the corresponding bremsstrahlung 
integrals. 



6.1 Resummation 

Before discussing the structure of the non-factorizable diagrams it is useful 
to recapitulate the resummation of the factorizable diagrams. The separation 
between resonant and non-resonant terms in the factorizable diagrams was made 
implicitly in section || by writing the corrections to production and decay as 
W(p 2 ) = W(m 2 ) + [W(p 2 ) — W(m 2 )], where only the propagator multiplying the 
on-shell contribution was resummed. This had the advantage that the second 
term is needed to one order less, as it is suppressed by a power of T ~ a. 
The first term is gauge invariant (this will be shown in section 0) . This scheme 
breaks down when W(m 2 ) contains logarithmic divergences (a charged unstable 
particle). The solution is to write the amplitude as a sum of fundamental 
loop integrals, the coefficients of which are gauge invariant. The logarithmic 
divergences occur in these integrals, so the resummation of these divergences 
does not influence the proof of gauge invariance. 

In the case of the non-factorizable diagrams for a single unstable particle, 
the resonant propagator now occurs inside a loop as l/({p + Q) 2 — m 2 ), with Q 
the integration momentum. Because of the infrared nature of the integral the 
main contribution is at Q ~ and we recover the original l/(p 2 — m 2 ), plus 
logarithms log(p 2 — m 2 ), up to the n-th. power for an n loop integral. These 
logarithms again render the residue at the pole p 2 — m 2 undefined, so once 
more we decompose the amplitude in a sum of fundamental loop integrals times 
coefficient functions. Now, however, the linear divergence is also contained in 
the fundamental integrals. A gauge invariant separation of the resonant terms 
is therefore accomplished by classifying the fundamental loop integrals into 
divergent and non-divergent ones for p 2 — > m 2 . 

The poles in the resonant fundamental loop integrals should of course resum 
to the same I/O 2 - M 2 ) as the propagator in the factorizable diagrams, with 
the arguments of the logarithms also shifting into the complex plane. Because 
of the logarithmic terms, we cannot perform the expansion (0) after integrating 
over the loops which comprise unstable particles. Just like in the case of the 



logarithmic on-shell divergences encountered in section 5A, the propagator can 
be resummed before this integration is done (in the fundamental integrals). 
Again, this expansion around the pole depends very much on the fact that 
this is the applicable kinematical regime; thus using the on-shell width is only 
valid when the integration momentum Q is not too large. As the singularity 
is infrared in nature, this holds for the first order (one loop); for higher orders 
more sophisticated techniques to isolate the poles will have to be used. 

In case of more than one unstable particle one again uses the same scheme for 
each propagator separately. The only difference is that now the non-factorizable 
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terms exhibit a much wider variety of singularities. The degree of divergence 
of the on-shell singularity obviously is equal to the number of unstable par- 
ticle propagators spanned, but the form is no longer the simple product of 
propagators that occurs in the factorizable terms. The one loop integrals are 
enumerated below. 

6.2 One loop integrals 

The first scalar function that can contain linear divergences is the four point 
function. Here we find three distinct singularity structures: one, two and three 
unstable particles in the loop. 4 Examples of these are shown in Fig. [?]. 




Figure 7: Examples of the three types of linearly divergent four point functions. 
The parentheses denote that the decay products are close to resonance. 

The first integral has been given in Ref. [?], also for a complex mass. The 
structure for the assignment in Fig. 01 is given by (A is the regulatory photon 



mass and —4§t = (t — m e 



D l = ^ 2 -- = 2 log 



in 



2\2 



Am^m 2 the Kallen function) 



s2 A 



t - (m A 



t - (m M - m e ) 2 + 2 V /3 ^ 



log 



'M 2 Z \ 



+finite terms 
(18) 



The introduction of the finite width does not pose any problems. 

The second four point function is more difficult to get correctly. The struc- 
ture is 



with (S1S3) 
by 



1 



/ <rSiS2S3S4 



(19) 



■p 2 + )/2, {s\Si) = {mly —p 2 _)/2 and the determinant given 
(sisi) ••• (s x s 4 ) 



1JS1S2S3S4 

u s 1 s 2 s a s 4 , 



(s 4 si) 



(s 4 s 4 ) 



(20) 



4 One needs one photon to obtain the singularity. There is one more divergent integral with 
no unstable particle propagator but regulated by a small photon mass A, which is unphysical: 
it diverges as 1/A (without logarithms). 
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(s iSj ) = (m? + m|-p|)/2 (21) 
(sip jJfe ) = {m}-ml-p%+pt k )/2 , (22) 

the momentum pj,- being defined as the difference of the momenta flowing 
through propagators with masses mi and rrij. The function /j are complicated 
functions of their arguments and the other on-shell invariants. 

The third type does not occur in W pair production and we have not yet 
investigated it. 




Figure 8: Example of a quadratically divergent five point function 



The five point function can maximally contain a quadratic on-shell diver- 
gence. If it is less singular, we just decompose it into five four point functions 
[?, ?] and apply the results given above. An example of a quadratically di- 
vergent five point function in W pair production is shown in Fig. In a 
quadratically divergent five point function the extra singularity occurs in the 
determinants multiplying the four point functions in the general expression 



En 



-1 



E w^rir a>(m + m + 2, » + 3) l (23) 



with the determinants 



SS1S2S3S4S5 



XSi Si + iSi+2Si+3 

u PiP2 Pa Pa 



(sisi) ••• (sis 5 ) 

(s 5 Sl) ••• (S5S5) 
(SiPl) ■■■ (SiPi) 

(si +3 Pl) ■ ■ ■ {Si+3Pi) 



(24) 



(25) 



Do(i,j, k, I) is the four point function with propagators with masses m,i,mj,mk 
and m;. We have to use complex masses in the overall factor, which is quadrat- 
ically divergent. This is reduced to a linear divergence by the coefficient de- 

in the four terms with a linearly divergent four point 



terminant 5 



■Si S i + lS i+ 2Si + 3 
P\P2 P3 Pi 



5 Again, in the physically uninteresting case that the all particles are taken on-shell and 
the divergence is regulated by a small photon mass the Eq diverges as 1/A 2 . 
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function. In the fifth term, where the photon is not included in the Dq, both 
the determinant and the four point function are regular, thus retaining the 
quadratic divergence. Note that, for consistency, we need a complex mass even 
in this last four point function, even if its main contribution is not at Q = 0.0 
This integral / d 4: Q/[(pi2Q)(pi3Q)(puQ)(pi5Q)] is regular, even though power 
counting would indicate a non-standard logarithmic divergence at Q = 0. 

All other higher point function are treated similarly. As we have had no im- 
mediate use for the cubically divergent six point function we have not analyzed 
it in detail. The quadratically divergent Fq is trivially decomposed in six five 
point functions. 

6.3 Soft bremsstrahlung 

As was noted above, the soft bremsstrahlung integrals {E~ < to <C T) have the 
same resonance structure as the virtual ones. The treatment given above can 
thus be taken over completely. It has been argued that the resonant part of the 
non-factorizable virtual diagrams (times lowest order) is equal and opposite to 
that of the corresponding soft bremsstrahlung graphs [?], so that the sum is 
non-resonant and can be neglected. This effect is also known from the Z pole 
[?]. However, as in the case of the logarithmic divergences, this cancellation 
only occurs when one integrates out the soft photons up to an energy w>T. 
Again, one cannot integrate when these small energy scales are relevant; of 
course if they are not important the narrow width approximation would have 
been sufficient. 



7 Gauge invariance 

In the previous sections we have resummed the factorizable resonant and non- 
factorizable resonant graphs in a complicated fashion. This gave rise to an 
amplitude of the form (for a single unstable particle) 

A =pT^M2+<P 2 )i (26) 

where M 2 , w and n(p 2 ) are computed from the analytic properties of the un- 
resummed amplitude. Now we show, essentially following M. Veltman [?], that 
this way of resummation does in fact give a gauge invariant result. Gauge 
invariance here denotes invariance with respect to any gauge group — the ar- 
guments are very general and only depend on the assumption that a fixed order 
amplitude is gauge invariant for all p 2 , which is used in the limit p 2 — ► m? 
(where the fixed-order calculation is not applicable for physics results). 

The first point to note is that the pole position M 2 is a property of the 
S matrix and a physically measurable quantity; it should therefore be gauge 
invariant. In this scheme it is given by the solution of the equation M 2 — m 2 — 
n(M 2 ) = 0. 

6 Unless a corresponding modification is made in the soft bremsstrahlung integrals. 
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The residue at the pole w is given as a sum of the resonant contributions 
from the factorizable and non-factorizable diagrams. However, the origins of 
the poles are quite different: in the non-factorizable diagrams it is provided 
by linearly (or higher) divergent fundamental loop integrals, in the factorizable 
diagrams by explicit propagators in the coefficients of these integrals. As the 
fundamental loop integrals are linearly independent, the resonant part of the 
non-factorizable diagrams, being composed of different fundamental integrals, 
must be separately gauge invariant to each order. 

For the factorizable part it is easiest to look at the unresummed amplitude 
truncated at a fixed order n, which is gauge invariant. In the (unphysical) limit 
p 2 — ► m 2 this diverges as l/(p 2 — m 2 ) n+1 ; the coefficient of this pole W- n -\ is 
therefore gauge invariant. Subtracting this pole and repeating the procedure 
for the remainder we find that all W-i are gauge invariant, in particular W-i, 
which gives the factorizable contribution to the residue in this order. This part 
therefore is also order- by-order gauge invariant. The occurrence of divergent 
logarithms does not spoil this argument due to the fact that they only occur 
inside the fundamental loop integrals. 

Finally, because the resonant parts of the unresummed amplitude are all 
gauge invariant, so is the sum of all the non-resonant contributions in the re- 
summed expression. 

8 Summary 

We have given a gauge invariant prescription to treat unstable particles in 
loop calculations, with explicit results for the one loop case. The scheme is 
an extension of the pole expansion method of M. Veltman [?, ?, ?] to include 
several charged resonances and loops. An advantage is also that only a minimal 
subset of diagrams is needed. It is shown that all necessary quantities can be 
obtained from a normal perturbative off-shell expansion. The problems which 
arise in the case of charged resonances with infrared divergences have all been 
addressed, as well as the corresponding treatment of the soft bremsstrahlung 
integrals (E^ < uj <C width). 

For a tree level calculation in this scheme one needs 

- the on-shell width U^\M 2 ) to one loop and 

- the resonant tree level amplitude W^\m 2 ). 
The 0{ol) corrections are given by: 

- improving the on-shell width II"( 2 )(M 2 ) to two loops, 

- the one loop factorizable diagrams W^im 2 ) with the substitution m 2 — > 
M 2 in the divergent scalar three point functions, 

- the resonant part of the one loop non-factorizable graphs, 

- the derivative terms W^°\m 2 )U(m 2 ) and W^'(m 2 )U(m 2 ), 



16 



- the tree level p 2 dependent resonant terms W^°\p 2 ) and 



- the tree level non-resonant terms. 



The soft bremsstrahlung integrals are treated analogously. 

The resonance structure of the non-factorizable graphs (those with unstable 
particles in the loops) is quite complicated in the case of multiple resonances 
(like W pair production): they are not just products of poles. These give a 
contribution which is not suppressed by additional powers of the width. The 
corresponding soft bremsstrahlung integrals only cancel this contribution if one 
integrates out the photons up to an energy much larger than the width. 

The expansion around the poles underlying our scheme leads to difficulties 
near (especially below) the thresholds for the production (or decay) of the un- 
stable particles. Connecting the regions above and below threshold will require 
other techniques. 



Acknowledgements. We would like to thank Fred Jegerlehner, Wolfgang 
Hollik, Thomas Sack, Andreas Salathe and Jos Vermaseren for many helpful 
discussions. One of us (GJvO) would like to thank FermiLab for the hospitality 
while debugging the complex functions. 



A Proof of Eqs (0) and (g) 

In order to prove Eq. (Rj) we write (writing s = p 2 ) 



n=0 



(27) 



We reexpress each term of the sum by expanding around s = M 2 and using 
m 2 — M 2 = -n(M 2 ) 



00 1 00 1 d nJrk 



ErrE^:L+fc lwm n (s)] s=M2 (-u(M 2 )) k 



n n\f-^k\ ds n+k 
n=0 fc=0 



E 

n,k=0 
oo 



1 (n + k\ d n+k r 
(n + k)\ I n I ds n+k 



2\\k 



W(s)U n (s){-U{M 2 )) 



=Af 2 



1 d ri 



n'=0 



2\\n' 



^^n.d^ W(s)(U(s)-U(M 2 )) 



s=M 2 



(28) 



with n! = n + k. Now we can expand n(s) as well as VF(s) around s = M 2 . 
Only the terms without any (s — M 2 )-factors will survive; therefore we get 



W{M 2 



£ w(M 2 )(n'(M 2 )f = ^ 



n'=0 



(M 2 ) 



w. 



(29) 
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The proof of the second relation (||) follows along similar lines. We start by 
writing N(s) as 

d n+e 



N(s) = £( 
1=1 

oo 

= E( 



s — m 



2\£-l 



E 



1 



(n + £)lds n+£ 



[W(s)H n (s)] £ 



n=0 
( oo £-1 



s — m 



E-E B 



1 <T 



n'! <is n ' 



n'(«) 



«=1 kn'=0 n'=0> 

with ra' = n + £. Using the previous proof the first term in braces gives 



E( 



s — m 



,£-iW(M 2 ) 



1 



W(M 2 ) 



1 



IF(M 2 ) 1 -n'(M 2 ) s — M 2 1 — n'(M 2 ) 



(30) 



(31) 



The second term is a simple Taylor series after changing the summation index 
tot' = £- re': 



E E ( s - »" 

t'=\ n'=0 



2\£'+n' 



oo oo 



E E( s -'"'' 

£'=ln'=0 



.! 1 d" 
re'! ds n ' 



W(a) 



2\r-i 



W(s) 



lF(s) s-m 2 -n(s) 
The sum is thus exactly n(s), which was to be shown. 



(32) 



B Frame dependence 

We will show that, up to terms of 0(ot 2 ), Eq. |3] gives a result which is indepen- 
dent of the definition of the other kinematical arguments of W. In particular, 
when these are angles, the result will only depend on the frame chosen through 
terms proportional to 0{a 2 ). 

Let us compute first in one frame, in which W is a function of p 2 and another 
set of variables Xj. The contribution of the factorizable diagrams than reads, 
up to 0(a), 



A 



(o+i) _ W(p 2 ,x i )-W(m 2 ,x i ) W(m 2 ,x l ) + ^ I W(p 2 ,x i )U(p 2 )\ p2 __ 



fact 



+ 



P 



M 2 



(33) 

In another frame we have Xi = Xi(yi,p 2 ), where the yi denote the kinematical 
variables in that frame, and the dependence of the Lorentz transformation on 
p 2 is shown explicitly. Suppressing the yi the amplitude now is 

^(o+l) _ W(p 2 , Xi (p 2 )) - W{m 2 ,x i {m 2 )) , W(m 2 , Xi (m 2 )) + ^ W(p 2 , x l (p 2 ))U(p 2 )\ p2=m2 



fact 



p2 _ m 2 



p 2 _ M 2 



(34) 



Substituting back Xi = Xi(p 2 ) the difference is given by 

, ai _ A \(0+l) _ W(m 2 ,Xj(p 2 )) - W(m 2 ,Xj(m 2 )) | W(m 2 , Xi (m 2 )) - W(m 2 ,Xj(p 2 )) 

V"' 4 A) fact — " + 



p2 _ m 2 



p 2 _ M 2 



+ 



P 



M 2 



.(35) 



18 



Expanding this around p 2 = m? and using M 2 = m 2 + n(m 2 ) + 0(a 2 ) this 
leaves only terms of order a 2 and (p 2 — m 2 )a. 

C Gauge invariance of the coefficients of scalar func- 
tions 

The fact that the coefficients of the scalar functions must be gauge invariant 
can be seen for the electroweak section of the standard model when we consider 
it with masses for all particles (except the photon and gluon). This does not 
break gauge symmetry. One can then consider the scalar functions in turn, and 
consider whether they can be linearly dependent. Note that this implies that 
the coefficients of scalar functions involving gauge-dependent parameters (for 
instance ghost masses) must vanish. 

- The one point function 

A (m 2 ) = -m 2 log(m 2 ) + ... (36) 
is obviously linearly independent of any other one point functions. 

- The two point functions 

' ^ ~ H l °< m 2 -m 2 + p 2 -^X ) + • • • 

(37) 

is unique for i < j and non-zero masses as the factor y/—6 does not occur in 
the one point functions or coefficients of other two point functions (except 
with permuted arguments, in which case the logarithms are different). 
The scalar two point functions with p\ = 0, m 2 = m 2 or mf = p 2 , m 2 = 
reduce to one point functions; the remaining possibility has a unique 
logarithm. 

- The three point function 

1 6 

C (m 2 ,m 2 ,m 2 k ,p 2 ,p 2 n ,pl) = — =^Li 2 (c<) (38) 

y-$(Pl,Pm,Pl) t=l 

must again be brought into a unique form i < j < k (if equal, order 
l,m,n). In case the kinematical determinant 5(pf ,p 2 n ,p 2 n ) is not equal to 
zero, it contains dilogs (or double logs), which are independent of any one 
or two point functions. The overall factor, again with a root, uniquely 
identifies the momenta; but different three point functions with the same 
external momenta but different masses will have different arguments for 
the dilogarithms. When the kinematical determinant is zero this no longer 
is the case, and three point functions do in fact reduce to two point 
functions [?]. Except for threshold (which we do not consider) these occur 
in counterterms (for instance jWW at q 2 = 0). 
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- The four point function depends on 6 masses and 10 momenta squared. 
It has the form 



with the Si internal momenta defined before. The overall factor has the 
same symmetries as the Dq and depends on all the masses and momenta 
squared in the box, so after bringing it to a standard order it uniquely 
identifies the four point function when it is non-zero. There are two 
cases where it is zero: real singularities on the boundaries of phase space 
{Sp\plpl = 0; this again occurs on threshold and in counterterms) and 
artificial ones within phase space. The latter ones depend on the values 
for s and t, at least one of which can normally be chosen freely. It can 
therefore not lead to a dependency which is valid for all allowed s and t. 

- The five point function is a linearly dependent set of five four point func- 
tions. However, this decomposition crucially depends on the dimension- 
ality of space (it uses the Schouten identity). As the gauge invariance 
of the amplitude does not depend on this, the coefficient of the Eq must 
be gauge invariant on its own. The same holds for six- and higher point 
functions. 

We have thus shown that all scalar functions with non-zero kinematical deter- 
minants are linearly independent. The others should be reduced to lower point 
functions. 




(39) 
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